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We develop an analytical theory of adhesive interaction of viscoelastic spheres in quasistatic approxima- 
tion. Deformations and deformation rates are assumed to be small, which allows for the application of 
the Hertz contact theory, modified to account for viscoelastic forces. The adhesion interactions are de- 
scribed by the Johnson, Kendall, and Roberts theory. Using the quasistatic approximation we derive the 
total force between the bodies which is not sufficiently described by the superposition of elastic, viscous 
and adhesive contributions, but instead an additional cross-term appears, which depends on the elastic, 
viscous and adhesive parameters of the material. Using the derived theory we estimate the contribution 
of adhesive forces to the normal coefficient of restitution and derive a criterion for the validity of the 
viscoelastic collision model. 
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1 INTRODUCTION 

Numerous phenomena observed in granular sys- 
tems, ranging from sand and powders to granular 
gases in planetary rings or protoplanetary discs, 
are direct consequences of the specific particle 
interactions. Besides elastic forces, common for 
molecular or atomic materials (solids, liquids, and 
gases), colliding granular particles exert also dissi- 
pative forces. These forces acting between contact- 
ing grains give rise to unusual properties of granu- 
lar matter. Hence, the use of an appropriate model 
of the dissipative interactions is necessary for the 
adequate description of granular systems. In real 
granular systems the particles may have a compli- 
cated non-spherical shape, differ in size, mass and 
material properties. In what follows, however, we 
assume that granular particles are smooth spheres 
of the same material. We also assume that particles 
interact exclusively via pairwise mechanical con- 
tact. 



There are three different types of forces act- 
ing between contacting particles: (i) repelling 
elastic forces, due to the compression of parti- 
cles, (ii) attractive adhesive forces which appear 
when particles share a common surface and (iii) 
dissipative forces, acting against the relative 
motion of the particles. The dependence on the 
material parameters and on the quantities which 
characterize the relative position and motion 
of particles is known for all of these forces, 
e.g. (Brillia ntov and Poschel 2()()4al) . The total 
force acting between the particles is, however, 
not just the sum of the above three compo- 
nents. Instead, a more careful analysis, sketched 
below, reveals that the total force contains an 
additional cross-term which depends on both 
dissipative and adhesive parameters. Using the 
obtained expression for the total force acting 
between adhesive viscoelastic spheres, we es- 
timate the effect of the adhesive force on the 



coefficient of restitution and analyze the range 
of validity of the frequently used viscoelas- 
tic interaction model dBrilliantov et al. 1 996 ; 
Brilliant ovand Poschel 2()()4bD . 

2 FORCES OF GRANULAR PARTICLES 
Elastic force. Consider a static contact of two 
spheres of radii R\ and R 2 . When the spheres are 
squeezed, the material in the bulk is deformed (Fig. 
[l]). The displacement field u (r ) causes the defor- 




Figure 1 . Collision of spheres, characterized by the time de- 
pendent compression = Ri + R2 — | r*i (t ) — r*2(i)| and 
its rate = «i(t) - v 2 (t). 



mation field 
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, hj = {x,y,z}. 



(1) 

In the elastic regime Uij(f) is proportional to the 
stress tensor a y '(f), which gives the z-component 
of the force, acting on a unit surface normal to the 
direction j (Land au andLifshitz 1 965): 



a l3 (r)=E 1 [u ij (r) - -5ijUn(r )\ + E 2 5 ij uu(r) . 

(2) 

Repeated indices imply summation and the elas- 
tic coefficients E\ and E 2 are related to the Young 
modulus Y and the Poisson ratio v by 
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The contact problem for elastic spheres has been 
solved by Hertz (118821) : the circular contact of ra- 
dius a gives rise to the elastic force, which read in 



terms of the compression £ = Ri + R 2 — \ fi — r 2 

a 2 = R e ff £ , 



where 
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The normal pressure P H {x,y) = a zz (x,y,z = 0), 
which acts between the compressed bodies in the 
plane of contact z = reads (Hertz 1882) 
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Adhesive force. The Hertz theory for the elas- 
tic contact of spheres was extended to adhesive 
contact by Johnson, Kendall, and Roberts (JKR) 
(1197 11 ). They found that the contact area is en- 
larged owing to the adhesive force and, thus, in- 
troduced an effective Hertz load Fh which would 
cause this enlarged area. The contact area of ra- 
dius a corresponds then to the compression £ H for 
the Hertz load F H . In reality, however, this contact 
radius occurs at the compression £ < The dif- 
ference between the Hertz compression and the 
actual one, £, was attributed to the additional stress 
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which is the solution of the classical Boussinesq 
problem (Timoshen ko 19701 ): This distribution of 
the normal surface traction gives rise to a constant 
displacement over a circular region of an elastic 
body. The displacement £ B corresponding to the 
contact radius a and the total load F B are related 
by 
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The quantity F B stands for the additional attractive 
adhesive force, which acts against F B . In the JKR 
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theory it is expressed by the adhesion coefficient 7, 
which is twice the surface free energy per unit area 
of a solid in vacuum: 



F B = 27ra' 



37 
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(9) 



Thus, the total force F = F H — F B and the actual 
compression, £ = £ H — £ B , read 
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implying a finite contact radius for F 

ajj = 67rD 7 i&. 
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(10) 



(11) 
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For decreasing force (F < 0) the contact radius de- 
creases until the minimal value 



sep 



(13) 



corresponding to maximal (in absolute value) neg- 
ative force. At this point the particles separate. 

Viscous force. If the deformation of contacting 
spheres changes with time, an additional dissi- 
pative force arises and the stress tensor con- 
tains an additional dissipative component o] 
For small deformation rate 
(ILandau and Lifshitz 1965 ), 



u xj \r , 



dis' 

it reads 



(*) = Vx 



(14) 

where 771 and 7] 2 are the viscous constants. 

If the impact velocity of the colliding bodies is 
much smaller than the speed of sound in the par- 
ticle material and if the characteristic relaxation 
time of the viscous processes in the bulk of the 
material is much smaller than the duration of the 
collision, one can apply the q uasistatic approx- 
imation dBrilliantov et al. 199 6). In this approxi- 
mation the displacement field u(f) coincides with 



that for the static case u e i(r), which is the solu- 
tion of the corresponding elastic problem. The field 
« e i(r ) in its turn is completely determined by the 
time-dependent compression £, i.e., u c \ = u e i(f, £) 
(Brilliant ov et al. 1996I) . Relating £ and a via Eq. 
(ITOb and neglecting a small hysteresis in the very 
beginning of the contact (Smit h et al. 19 89). we 
obtain for adhesive interaction 

• d d 

u(r,t) ~ £^TM e i(r,£) = a-g^Uei (r,£(a)) . (15) 

The dissipative stress tensor reads, respectively 
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From Eqs. (fToT) and © follows the relation be- 
tween the elastic and dissipative stress tensors in 
quasistatic approximation (Brilliantov 2005), 



dis 



Vl, E 2 



V2) 



(17) 



meaning that the dissipative tensor is obtained 
from the corresponding elastic tensor by substitut- 
ing the elastic constants by the viscous constants 
and applying the operator ad /da. In particular the 
normal component of the stress tensor at the plane 
z = for an adhesive contact reads 

a*?(x,y,z = 0) = P H (x,y) - P B (x,y) , (18) 

with P H (x,y) and P B (x,y) given by Eqs. (161719b . 

From Eqs. (1 1 VII 8b we find the dissipative stress 
at the contact plane and, integrating it over the con- 
tact area, the dissipative force. Referring for detail 
to ([Brilliantov 2005 ), we present here the final re- 
sult: 
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The first term in the rhs of Eq. dT9l corresponds 
to the dissipative force in the absence of adhesion. 
The second term is the corresponding cross-term, 
which depends on both adhesive and dissipative 
parameters. 

Remind that the above relations for the dis- 
sipative force have been obtained within the 
simple JKR theory, which is quite accurate in 
the range of param eters of practical interest 
dAttard and Parker 1992b . 

3 COEFFICIENT OF RESTITUTION 

An important characteristic of rarefied systems 
(granular gases (Brilliantov and Posche l 20041)10 
is the coefficient of restitution, which quantifies 
the loss of mechanical energy for pairwise colli- 
sions. It relates the pre-collision relative velocity, 
g = v\ — f 2, to that after the collision, g' = v[ — v' 2 : 



e(g) = -g'/g- 



(22) 



The coefficient of restitution may be evaluated 
solving the two-body collision problem with given 
interaction forces, yielding g' as a function of g. 

For small velocities, when the kinetic energy of 
the relative motion of colliding particles is close 
to the surface energy of the contact, the adhesive 
forces may change the coefficient of restitution 
qualitatively. Indeed, adhesive particles may stay 
compressed in contact even if the external load 
vanishes. That is, a tensile force must be applied 
to separate the particles. The work against this ten- 
sile force at the very end of the collision reduces 
the kinetic energy of the relative motion after the 
impact, that is, it reduces the coefficient of restitu- 
tion. For small impact velocity the kinetic energy 
of the relative motion may be too small to over- 
come the attractive barrier, i.e., the particles stick 
together after the collision corresponding to e = 0. 

Consider first a pure viscoelastic collision with 
the impact velocity g and coefficient of restitution 
e v . Using the definition (|22b . the energy balance 
reads 



\m^g 2 - \m^(g)g 2 



W, 



dis • 



(23) 



where 



eff m x m 2 
m = 



(24) 

mi + 777,2 

The work W& s results from the dissipative 
force Eq. (fT9l with 7 = 0. The correspond- 
ing coefficient of restitution e v is analyti- 
cally known ( Schwager and Poschel 1998 ; 
IRamfrez et al. 1999b . Turn now to viscoelastic 
impacts with adhesion, characterized by the 
coefficient of restitution e ad . The energy balance 
reads 

I m ey _ 1 rrf^gy = WdisA + w disB + Wad , 

(25) 

where WdisA is the work of the dissipative force 
F dis due to the first term in Eq. (fT9l) and WdisB 
is the work due to the second term. Finally, W ad 
is the work due to adhesion, i.e., it results from 
the adhesive force in the region where the total 
force Eq. (fTTT) is negative, that is, in the region 
where the contact radius varies from ao to a sep 
(Brilliantov and Poschel 2004al) : 
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Using the approximation W d [ sA ~ from Eqs. 
(123125b we obtain the coefficient of restitution for 
viscoelastic collisions with adhesion: 



4i(£) = 4(9)- 



2(W ad + W disB ) 

m efig2 



(27) 



W ad may be found employing Eq. (fTTT) for the total 
force F(a), Eq. (ITOt for the compression, which 
allows to obtain d£ / da, and Eqs. (1 1 21 1 31) for a and 



l sep- 



5„,5 n 2 d4 \ 1 / 3 



W ad = q (n 5 1 5 D 2 Rt s ) 



(28) 



with the analytically known pure number q ~ 
0.578 (IBrilliantov and Poschel 2()()4al) . For small 
dissipation and small adhe sion, Wdis B may be 
roughly estimated as (Brillia ntov 20051) : 



disB 



qxB^D-Wg*'* (m eff R 2 eff ) 3/10 



(29) 
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where q x m 7.93. 

From Eq. d27t we obtain the condition for the 
validity of the viscoelastic collision model: 

# 2> 4ff (Wad + WdisB) , (30) 

which with Eqs. (128129b may be written as 

9 2 »9 2 c = (l 5 D 2 Rt s ) 1/3 /m^ (31) 

B « (D 2 hR 2 en f 6 (m eff ) 1/2 . (32) 

With the above condition for B one can neglect 
W-disB as compared to W a d for the impact veloc- 
ity g being of the order of g c {g ~ g c ). Then we 
obtain respectively the condition for sticking colli- 
sion when £ a d(fi' s t) = 0: 

^m ef! e 2 v (g st )gl = W ad . (33) 

Hence, if g < g st for head-on collisions (vanishing 
tangential component of the impact velocity), the 
colliding particles stick together and form a joint 
particle of mass mi + m 2 . 

4 CONCLUSION 

The collision of adhesive viscoelastic spheres is 
characterized by (i) the elastic Hertz force, (ii) the 
dissipative force originating from viscoelastic bulk 
deformation, and (iii) the adhesive force. We use 
the continuum model of adhesive contact by John- 
son, Kendall, and Roberts (119711) which is ade- 
quate in the range of parameters of practical inter- 
est. The total force was derived under the approx- 
imation of quasistatic deformation, that is, the im- 
pact velocity is assumed to be much smaller than 
the speed of sound in the material and the viscous 
relaxation time is much smaller than the duration 
of the collision. This force is not only the superpo- 
sition of its three components (i-iii), but there ap- 
pears an additional cross -term, which depends on 
both viscous and adhesive parameters of the mate- 
rial. 



Using this force we estimated the contribution of 
adhesive forces to the normal coefficient of resti- 
tution as well as the range of validity of the vis- 
coelastic collision model (Brilliantov et al. 1996) 
and the condition for sticking impact of head-on 
collisions. 
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